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I— I i Abstract 
^ ■ 

Qh| We investigate the interplay between crossed channel final state interactions and the con- 

straints from two-particle unitarity for the reactions J/iJj — > Virn and VKK, where V is either 
. Lo or (/). Using a model where the parameters are largely constrained by other sources, we 

find that, although small, crossed channel final state interaction can influence the amplitudes 
^ ' considerably, in special areas of phase space. These results cast doubt on the inapplicability of 

^ : 

QQ , unitarity constraints on production amplitudes as recently claimed in the literature. 
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1 Introduction 



In recent years, much interest and effort have been put into extracting information about two- 
body interactions from their Final State Interaction (FBI) in production reactions. To study the 
. final state interaction among multiple final particles rigorously, one in principle needs to solve 
complicated many-body equations, where not all of the necessary input is known. However, if one 
of the pair-wise interactions in the final state is much stronger than all the others, one may use a 
simplified method, first introduced by Watson [Ij. The theorem may be derived from the unitarity 
relation for two-particle interactions from a production reaction 

lm{A,) = -akT*„Ak, (1) 

where Af^. is the full production amplitude into channel k, is the two-body phase space factor for 
channel k and T*^. is the complex conjugate of the scattering amplitude connecting channels i and 
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k. Here summation over k is assumed. This relation is also called Extended Unitarity (EU) [2j and 

has been used in many analyses. We will also use this phrase in what follows. The case studied by 

Watson was the single channel case. Since Im(j4) is a real number, it automatically follows then 

that the phase of the production amplitude equals to that of elastic scattering — this is known 

as the Watson theorem. Then the final state interaction is well described by the Jost-function for 

the scattering amplitude [3], which for large scattering lengths and small momenta agrees with 

the scattering amplitude. However, in many applications, as the J/'0 decays we will study in this 

work, the mentioned conditions are not met in the whole kinematic regime. For example above 

the KK threshold coupled channels become relevant. Since Eq. ([T|) contains a summation, there is 

no connection anymore between the phase of the scattering amplitudes and that of the production 

amplitudes. In addition, there might be effects from crossed channels. Those were most recently 

discussed in Ref. [1]. The effects of coupled channels on the mentioned J/^p decays were first 

investigated in Ref. [5]. The idea was refined in Ref. [6j and later improved in Ref. [7] — especially 

here the data of Refs. [SlEllin] was included. Crossed channel effects were studied in Refs. [TT|[T^. 

however, no detailed comparison to the most recent data was performed. 

Recently, Bugg [13J argued that the EU is in confiict with experimental data when applied to 

the case with multiple resonances being present in one partial wave, such as the scalar meson 

channel. Based on the incorrect assumptiorj^ that Watsons theory was originally derived without 

overlapping resonances, the author argues that the extension of this theory to describe the scalar 

meson production, in which /o(600) (the so-called a meson) and /o(980) are present, is not reliable. 

The central statement of Ref. [13] is that, since EU predicts the same phase of production amplitude 

and scattering amplitude (and he assumes this to be the case even above the opening of inelastic 

channels), analyticity would further require that the relative magnitude of the two resonances 

should also be the same for both processes. He then demonstrated that this is in conflict with the 

PWA results of the s-wave contribution to the vrvr invariant mass spectrum in the J/ip ^ wvrvr, 

where there is no signal of a deep dip around the KK threshold as it would follow from the reasoning 

presented above. In this work we investigate, if this observation indeed shows that Eq. ([T]) is wrong. 

Before addressing this issue, one should notice that the naive EU as given in Eq. ([T|) only tells us 

the relation between the production amplitude which only includes the information of ESI of two 

final particles and the corresponding elastic scattering amplitude of the two particles. This means 

^The reasoning by Watson is given for a single continuum channel. No statement is made on the number of 
resonances allowed. 
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that the FSI between these two particles and other hadrons in the final state, which will be called 
crossed channel FSI in the following, has been neglected. Sometimes, such an approximation can 
be employed. For instance, for the dipion transition between charmonia tp' J/ipnir, if there is 
no resonance which couples strongly to J/ipir or ip'n. In addition, the interaction between J/^lJ and 
ijj' and the pions is OZI suppressed and therefore small [TJj. However, the situation is different 
for the J/V' decays into three light flavored mesons. For the decay J/ip ^ lotttt, there is a well- 
established axial vector meson, the 5i(1235), which couples strongly to both the ojtt and J /ipir [15j 
(certainly the coupling strength of biwvr is much larger than that of J/tpbin since the latter is OZI 
suppressed). 

It should also be stressed that a linear relation like Eq. ([T]) is putting a lot weaker constraints 
on the amplitudes as a non-linear relation like the corresponding one for elastic scattering does. 
In particular, the amplitude gets fixed only up to a polynomial. Thus, it is not correct that 
the requirement of unitarity and analyticity can determine the relative strengths from different 
resonances as well. On the contrary, based on Eq. ([T]) there is still the freedom to adjust the 
relative coupling between different resonances. 

The main aim of this work is to demonstrate that the conflict between the experimental data 
analysis and the result that follows from EU can be largely removed by the crossed channel effects. 
Here we mainly concentrate on the production of one vector meson and two pseudoscalar mesons in 
J ftp decays where the two pseudoscalars are in relative s~wave. We follow the ideas and formalisms 
in Refs. [6l[7l[ll]. In Ref. [7], the authors constructed the production amplitude from two pieces. 
The first part is the coupling of the J/^p to a vector meson and a scalar source, and the second 
part is the coupling of the scalar source to two pseudoscalar mesons. The latter is described by 
scalar form factors. Final state interaction is incorporated through a chiral unitary description 
with coupled channels, in which the scalar mesons (cr and /o(980)) are produced dynamically from 
the s-wave interactions between the Goldstone bosons. In this way, it is possible to investigate 
the OZI violating effect and to extract the low energy constants (LECs) of chiral perturbation 
theory (ChPT), if scalar form factors are calculated up to next-to-leading order (NLO). In Ref. [7] 
a deep dip appears right above the KK threshold in the s-wave vrvr invariant mass spectrum of the 
J/ip ^ iOTTTT as discussed above. However, this is apparently in conflict with the PWA results by 
the BES Collaboration [S]. Bugg considered this as a signal that EU is wrong. However, it should 
be noted that in Ref. [7J , the OZI violation parameter is assumed to be real. This means that the 
contribution of crossed channels was neglected. In fact, in Ref. [TT] the authors have investigated 
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the contribution from those mechanisms. In their work, the crossed channel FSI are modelled by 
some three-meson loop diagrams. It is shown that the contribution of those mechanisms is not 
negligible. It is clear that the crossed channel FSI can distort the prediction of the naive EU. After 
considering those diagrams, it is possible to largely remove the conflict between the prediction of 
EU and the BES PWA results as will be shown in this paper. 

The formalism of our analysis in presented in section [2j In section [3l we give the results of our 
analysis and the parameters from fitting to the data. The conclusion and some discussions are 
given in sectional 

2 Formalism 

The experimental data of the J/V' decaying into a vector meson and two pseudoscalar mesons 
are collected in a series of papers [SI El [13 [13 [IZ] • The latest data are published by the BES 
Collaboration in Refs. [H [Qj [TO] where also a partial wave analysis (PWA) is provided. In this 
way the s~wave contribution of the final tttt and KK pairs was isolated from the other partial 
waves. It is convenient for us to use these PWA results because we can concentrate on the scalar 
meson channel and avoid the complication from the contribution of other resonances. Hence, our 
analysis is based on the PWA results [HI [9]. However, it should be stressed that vrvr s-wave in 
the region of interest is only a very small part of the full signal which emerges from interferences 
of this partial waves with others. As a result, the extracted s-wave strength should depend on 
its assumed phase motion. We expect that the amplitudes of our analysis will show a different 
phase motion compared to those used in the analysis so far. In addition, in the BES analysis some 
possible tree level resonance contributions were neglected, such as the subthreshold contributions 
from the p (for wttvt) and the K* (for (f)KK), and the contribution from Ki mesons (for (j)KK) as 
shown in Fig. [TJ It has been stressed in Ref. [18] that the contribution from the p is important and 
may have significant influence on PWA result. This uncertainty also enters our work because our 
analysis is based on BES PWA result. Because of these reasons we do not aim at a high quality 
description of the data, but only at pointing at possible deficits of previous analyses, especially of 
Ref. [13j. Eventually the full data set should be reanalyzed based on improved amplitudes. 

Usually in PWA only some tree level amplitudes are constructed, and loop contributions are 
effectively absorbed into the coupling constants. This is the reason why in some phenomenological 
analyses the coupling constants are allowed to take complex values. But it should be noted that in 
Ref. ^ the OZI violation parameter and coupling constant are set to be real. This implicitly means 
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Figure 1: The tree level Feynman diagrams included in our work for the decays J ^ wvrvr and 

that the crossed-channel contributions are not included. However, the diagrams shown in Fig. [2] 
can contribute as stated in the Introduction, if the outgoing light meson pairs are in the s-wave. 
The formalism we will use to evaluate the mentioned mechanisms follows Refs. [71[TT]. Here we 
just illustrate the formalism briefly, and the detailed description can be found in Refs. [7} lllj. 

The Feynman diagram that only includes the FSI between the pseudoscalar mesons is given in 
Fig. El Using the definition and normalization of Ref. the corresponding amplitudes can be 
written as 



-M^" = y|c<^(0|ss + A^nn|7r7r; 



7=0 



)^C^{^\s-s^\^nn\KK) 



1=0 



(2) 
(3) 



Here Ai'^'^ and Ai^^ are the amplitudes for the processes J/t/j ^ ipiTir and J/t/j ^ (pKK, respec- 
tively. Xfjj is the OZI violation parameter. The amplitudes for J ftp ajnir and J/ip ^ uoKK can 
be similarly written as 
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with 



A., 



C^{0\ss + X^nn\KKy='^. 
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(5) 
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(7) 



The scalar form factors are defined as 
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{0\nn\KK)^=°. 
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J/HJ K'(K-) W J/tp K,(K^) W 

Figure 2: Three-meson loop diagrams for J/ip ^ ivmr. (The corresponding diagrams for other 
channels can be obtained by changing to to (f) and vrvr to KK in the last two diagrams(for J/-^ — > 
(f)KK) or changing u; to i?!> in the last two diagrams(for J jil) ^ ^vrvr) respectively.) 



Here Bq denotes the strength of the scalar quark condensate. The scalar form factors can be 
calculated within ChPT to a given order. However, because we are interested in the energy 
range up to 1.2 GeV, where the energy is too high for ChPT to be applied, the Chiral Unitary 
Approach IS] is adopted. This method is described in detail in Ref. [19J, and the expressions 
for the scalar form factors up to NLO can be found in Ref. [7]. 

As discussed in the Introduction, crossed channel effects are potentially important. We model 
them by some three-meson loop diagrams following Ref. [1 IJ . Those diagrams are shown in Fig. [5J 
The effective Lagrangians for the corresponding vertices can be written as [TT] . 

CAiB)VP = D{V^AB^\P})-iF{V^,[A^\P]) 

Ca(b)^p = Dij^uiB^'P) (10) 

where (• • •) means SU(3) flavor trace and V{P) are the SU(3) matrix representations of the vector 
(pseudoscalar) mesons as usual. Using the tensor field formalism to describe the spin-1 mesons |20j . 
A^^'^ and B'^'^ are the SU(3) matrices for axial- vector mesons with J^*^ = and 1"' , respectively, 
and V^'^ is the SU(3) matrix for vector mesons. 
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Figure 3: Feynman diagram of the J /ij: decays into one vector meson and two pseudoscalar mesons 
with the FSI of two pseudoscalar mesons. 



The propagators of vector and axial vector mesons in tensor formalism are given by [20] 



iV 



fiupcr 



= {0\T{W^,WpM 



' fiiy " pa i 



9pp 9ua {M^ - P^) + 5mp Pa - 9pa Pv Pp-if^^ z^)l , (11) 



where Mw and P are the mass and momentum of the spin-1 meson, respectively. 

In the standard way, it is easy to write down the amplitudes of those diagrams. The squares in 
the diagrams mean the full meson-meson scattering amplitudes involving the loop resummation 
through the BS equation [19]. The three-meson loop integrals are divergent. They are regularized 
using a cut-off as described in detail in Ref. [21]. 

Using the amplitudes described above, it is straightforward to calculate the invariant mass 
spectra and fit the free parameters in the amplitudes to the date 



3 Numerical Results 

In our work, we will only consider the J/ip decays into the wvrvr, (j)7nT and (j)KK channels. Although 

it will be straightforward to include the ujKK channel as well, as pointed out in Ref. [3, we do 

not include this channel in the fitting procedure due to its large uncertainty. For illustration our 

results for that channel are shown below. The free parameters basically come from the diagrams 

shown in Fig. [3l which include the NLO LECs L4, L5, Lg and Lg appearing in the scalar form 

^Note, the data as given in the publications are to be corrected for an invariant mass dependent acceptance. See 
Ref. [2] for details. 
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factors, the parameter C,^ and the OZI violation parameter A,^. Other parameters such as C^^ and 
can be related to and through Eqs. |6|7| . The cut-off in the two-meson loop is set to 
be 0.9 GeV as Ref. |7j. Most of the parameters appear in the three-meson loop diagrams such 
as the coupling constants, the mixing angle between Ki mesons are fixed in Ref. [TT] and some 
related papers [52]. In order to reduce the number of free parameters, we adopt those values and 
keep them fixed in our calculation. The uncertainties due to those coupling constants are included 
in the error band in our results shown in Figs. [5]and[6l We adopt A = 1.0 GeV for the cut-off 
parameter used in three-meson loop calculation. Note that the quality of the fit does not change 
considerably when the cut-off is varied around that value — we varied A in the range between 0.85 
and 1.15 GeV; the resulting variation in the observables is smaller than the error bands shown, if 
the parameters are refit for each cut-off. Besides the magnitudes of the coupling constants, their 
relative signs are also important. Most of the relative signs have been fixed in Ref. [11]. But there 
is still a relative sign between the J/ipVP and J/ipAP couplings which cannot be determined. 
This leads to two sets of solutions in Ref. [TT], and thus we keep it as an adjustable freedom in 
our work. 

As mentioned above, we assume that the tree level contributions shown in Fig. [1] are effectively 
absorbed into the partial wave amplitudes of the BES analysis. So in our calculation, we include 
the contributions from p, K* and Ki exchange. 

Another difference between our work and Ref. [7] is that we use equal weights for different 
channels and fit the s-wave contribution up to 1.2 GeV for every decay including the ojtttt channel, 
while in Ref. [7] the fitting range stopped before the deep dip appearing in their description of the 
cjvrvr channel. 

With the considerations described above, we totally have 6 free parameters and a relative sign 
between the coupling J/ij)VP and J/ipAP, which we fit to the invariant mass spectra using the 
MINUIT program [23j. The results for the best fit and a comparison with the BES data are 
shown in Figs. [5] and El where also the uKK channel is shown for illustration, although it was 
not included in the fit. The resulting parameters are given in Table [TJ The uncertainties from the 
coupling constants and the free parameters are included in the error band. 

4 Discussion and Conclusion 

Before discussing the results in detail we illustrate the effect of the three-meson loop diagrams in 
some more detail. In Fig. [Hwe show a comparison of magnitude and phase motion for the two-kaon 
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Figure 4: Comparison of magnitude and phase motion of the two-kaon loop (dashed line) and the 
three-meson loop with an additional K* exchange (solid line). 

Table 1: Parameter values obtained by fitting to BES results for the J/ip decays into the lvittt, 
(j)7nT and (j)KK. All values of the renormalized LECs L\ are quoted at the scale /i = rup, and the 
relative sign between the J/ipVP and J/tpAP couplings is even (x^/d.o.f. ~ 2.1). The results of 
Ref. [7j is also given for comparison. 
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0.9±0.025 



loop (Fig. [3j) and a corresponding three-meson loop with additional K* exchange (third diagram 
of Fig. where for illustration the scattering T-matrix is replaced by a constant. As one can 
see, the three-meson loop leads to both a non-trivial variation in magnitude and phase compared 
to what emerges from the kaon loop alone. Especially, there is already a phase motion present 
for TTTT invariant masses below the two-kaon threshold. This is a result of the consideration of the 
width of the K* by which the kinematically allowed KKtt can contribute. Other three-meson 
loops show different patterns. Therefore crossed channel effects lead to contributions that cannot 
be included by just using complex coupling constants. 

Figures [5] and [6] show that we describe the data quite well, although the description near the 
KK threshold in lutttt channel where the deep dip originally appears is not good enough, it has 
been improved a lot compared to Ref. \X\ ■ The LECs we get are consistent with other works within 
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Figure 5: Solid line: full results for J/ip ^ Vtttt corresponding to parameters shown in Table [2 
dotted line: tree level contributions from p; dashed line: results without the contributions from 
the three-meson loop diagrams and the tree level diagrams. The BES PWA results for the s-wave 
contribution are shown by the histograms. 

uncertainties. However, the parameter which describes the coupling of J/V' to vector mesons 
and a scalar source is about a factor of 1.5 larger than that in Ref. [7]. This is needed to cancel 
some contribution of new mechanisms we added, otherwise we cannot fit the data well. 

In our work, we investigated the interplay between crossed channel effects and the prediction of 
EU. It seems from our calculation that the prediction of EU can be influenced by crossed channel 
FSI greatly in those regions of phase space where the amplitudes from the two-body interactions 
are small. So any conclusion about the conflict of EU with experimental data, which is based 
on those regions, should be questioned, because there might be other mechanisms besides those 
governed by EU that are relevant. Note, however, that overall the description from just keeping 
the two-body vrvr and KK interaction is quite impressive. To separate the influence from crossed 
channel effects from the leading amplitudes is difficult and may be model dependent. In our 
work, we follow the ideas of Ref. [TT] to calculate the crossed channel FSI. The calculation is 
model dependent and not complete in principle, but it offers some information about the interplay 
between EU prediction and crossed channel FSI. It seems that we cannot confirm a conflict of EU 
predictions and the experimental analysis. 

We also addressed uncertainties in the PWA by the BES collaboration. In their analysis some 
possible tree level contributions were neglected, which contribute in principle and are not negligible 
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Figure 6: Solid line: full results for J /ip VKK corresponding to parameters shown in Table [U 
dotted line: tree level contributions K* and Ki mesons; dashed line: results without the contribu- 
tions from the three-meson loop diagrams and the tree level diagrams. The BES PWA results for 
the s-wave contribution are shown by the histograms. 

in our calculation (see dotted lines in Figs. and With these uncertainties in mind, we conclude 
that there is no conflict between EU and experimental data. Especially, there is no justification 
to replace Eq. ([T|) by something else that has no theoretical ground, as was done in Ref. [13]. 
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